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ABSTRACT 

We provide a new estimate of the local supermassive black hole mass function 
using (i) the empirical relation between supermassive black hole mass and the Sersic 
index of the host spheroidal stellar system and (ii) the measured (spheroid) Sersic 
indices drawn from 10k galaxies in the Millennium Galaxy Catalogue. The obser- 
vational simplicity of our approach, and the direct measurements of the black hole 
predictor quantity, i.e. the Sersic index, for both elliptical galaxies and the bulges of 
disc galaxies makes it straightforward to estimate accurate black hole masses in carly- 
and late-type galaxies alike. We have parameterised the supermassive black hole mass 
function with a Schechter function and find, at the low-mass end, a logarithmic slope 
(1 + a) of ~ 0.7 for the full galaxy sample and ~1.0 for the early-type galaxy sample. 
Considering spheroidal stellar systems brighter than Mb = — 18 mag, and integrating 
down to black hole masses of 10 6 M©, we find that the local mass density of super- 
massive black holes in early-type galaxies Pbh, early— typo 

= (3.5 ± 1.2) x 10 5 /i? M o 

Mpc -3 , and in late-type galaxies Pbh.iatc-typc = (1-0 ± 0.5) x lO 5 /i 3 o M0 Mpc -3 . 
The uncertainties are derived from Monte Carlo simulations which include uncer- 
tainties in the M^-n relation, the catalogue of Sersic indices, the galaxy weights and 
Malmquist bias. The combined, cosmological, supermassive black hole mass density is 
thus fibh, total = (3.2 ± 1.2) x 10~ 6 hjQ. That is, using a new and independent method, 
we conclude that (0.007 ± 0.003)/i 3 per cent of the universe's baryons are presently 
locked up in supermassive black holes at the centres of galaxies. 

Key words: black hole physics — galaxies: bulges — galaxies: fundamental param- 
eters — galaxies: luminosity function, mass function — galaxies: structure — surveys 



1 INTRODUCTION 

Two purely photometric properties of galaxies, or rather 
their spheroidal components, are known to correlate 
strongly with a galaxy's supermassive black hole (SMBH) 
mass Mbh- The first property is optical luminosity (Kor- 
mendy 1993; Franceschini et al. 1998; Magorrian et al. 1998). 
Due to the observation that SMBHs are associated with the 
'bulge' of a galaxy, and not the disc, it is necessary to per- 
form a bulge/disc decomposition if one is to properly treat 
lenticular and late-type galaxies. At present, apart from (Er- 
win, Graham & Caon 2002; their figure 3, with only eight 
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1 By the term 'spheroidal', we mean an entire elliptical galaxy or 

the dynamically hot component of a disc galaxy. 



elliptical and five disc galaxies) , no opticallj0 calibrated rela- 
tion that pertains to both elliptical galaxies and the bulges 
of disc galaxies is available. On the other hand, one can 
simply exclude the disc+bulge galaxies and only work with 
elliptical galaxies (e.g., McLure & Dunlop 2002, 2004). 

In an effort to include disc galaxies, some authors have 
assigned some fixed fraction, such as three-tenths, of a disc 
galaxy's total light to that of the bulge. However, given 
there is a known trend of decreasing bulge-to-total lumi- 
nosity ratio with increasing morphological type (e.g., Hub- 
ble 1926, 1936; Kent 1985; Simien & de Vaucouleurs 1986; 
Andredakis, Peletier & Balcells 1995), the above approach 
introduces a systematic bias such that the SMBH masses 



2 A useful near-infrared M^h - ^spheroid relation is presented by 
Marconi & Hunt (2003), updated in Graham (2007). 
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are over-estimated in the late-type disc galaxies and under- 
estimated in the early-type disc galaxies, skewing the mid- 
to low-mass end of the SMBH mass function. Another ap- 
proach has been to use the average bulge-to-total flux ratios 
derived from past i? 1,/4 -bulge + exponential-disc decomposi- 
tions (Simien & de Vaucouleurs 1986). However, Andredakis 
& Sanders (1994) showed that Sb and Sc galaxies are, on av- 
erage, better described with an exponential-bulge than an 
R 1 ^ -bulge. Andredakis et al. (1995) subsequently showed 
that an _R 1//n -bulge was more appropriate, with the Sersic 
index n shown to decrease with increasing disc galaxy type. 
While an 7? 1 ^ 4 -bulge plus exponential-disc results in an over- 
estimate of the B/D ratio when the bulge has a Sersic index 
n < 4 (c.f. Figures 15 and 1 in Graham 2001), the use of 
an exponential-bulge + exponential-disc decomposition re- 
sults in an under-estimate of the bulge-to-disc (B/D) flux 
ratio when the bulge actually possesses a Sersic profile with 
n > 1 (Graham 2001, his figure 13). Therefore, 7? 1/n -bulge 
plus exponential-disc fits are required. The various M^-L 
relations in the literature predict SMBH masses that differ 
by factors of two to ten depending on the luminosity. This 
obviously inhibits the use of the Mhh-L relation at present. 
An investigation of this problem is not undertaken here but 
presented in Graham (2007). 

The second^ photometric quantity known to correlate 
with Mbh is the concentration of the stars in the host 
spheroidal stellar system (Graham et al. 2001, 2003). This 
concentration is monotonically related to the shape, i.e. the 
Sersic index n, of the spheroid's light-profile (Trujillo, Gra- 
ham & Caon 2001, their equation 6). Moreover, the M^-n 
relation is known to be as tightly correlated as the M^-a 
relation and have the same small degree of scatter (see No- 
vak, Faber & Dekel 2005 for a recent comparison of these 
relations) . 

Using an expanded galaxy set (27 galaxies) with up- 
dated distances and black hole masses, Graham & Driver 
(2007a) have recently shown that the M^-n relation (see 
FigO is curved rather than linear. Fitting a quadratic equa- 
tion, they obtained 

log(M bh ) = 7.98(±0.09) + 3.70(±0.46) log(ra/3) 

-3.10(±0.84)[log(n/3)] 2 , (1) 

with an intrinsic scatter extrinsic = 0.18to'oei dex. The to- 
tal absolute scatter in log M b h is 0.31 dex, which compares 
favourably with the value of 0.34 dex from the log Mbh-c 
data and relation in Tremaine et al. (2002). The parameter 
in front of the second order term in Eq[T]is inconsistent with 
a value of zero at the 99.99 per cent confidence level. 

We have explored here whether the departure from a 
linear relation may have been driven by an increased and 
uneven scatter, i.e. outlier^], at the high-mass end of the 
Mbh-n relation, or whether the curvature is inherent in the 
rest of the data set. In Figf2]we show the results of fitting 



3 A third photometric quantity that has been predicted to cor- 
relate well with SMBH mass is the central stellar density of the 
host spheroid (Graham & Driver 2007a, their section 6). 

4 As noted in Graham & Driver (2007a, their Section 3.3), there 
is reason to suspect that the highest SMBH mass, pertaining to 
NGC 4486, may have been overestimated, perhaps by a factor of 
4. 



a log-quadratic relation after the removal of the five highest 
mass data points from Fig[T] The coefficient in front of the 
quadratic term is again found to be inconsistent with a value 
of zero, this time at the 99 per cent level, and all three 
coefficients remain consistent, at the la level, with the values 
given in Eq[T] 

Support for the Mbh-n relation stems from its appli- 
cation to galaxies not included in its construction. For ex- 
ample, De Francesco et al. (2007) have recently measured a 
Sersic index n — 4.1 for NGC 3998, from which one would 
predict log(Mbh/MQ) = 8.43, in perfect agreement with the 
mass they derived from a kinematical study of the nuclear 
gas. In another example, Guhathakurta et al. (2006, their 
Figure S2) report a Sersic index n = 2 for M31, which has a 
SMBH mass equal to 3.5xlO 7 M (Ferrarese & Ford 2005) 
and is therefore also in good agreement with the data in 
Fig 03 Furthermore, for a sample of 11 narrow-line Seyfert 
galaxies (Ryan et al. 2007), the Mbh-n- relation predicts 
SMBH masses in agreement with those derived using the size 
of the broad line region and the continuum flux, and suggests 
a problem with the (Mbh-luminosity)-derived masses. 

The Mbh-n- relation also implies a maximum mass to 
which SMBHs have formed. The broad turn-over seen in 
FigH] peaks at n = 11.9, where the predicted la range of 
SMBH masses spans 0.8 to 3.8xlO 9 M (Graham & Driver 
2007a, their Eq.8). In Graham & Driver (2005, their fig- 
ure 1), one can see that for n >= 5 there is not that much 
difference in profile shape, and hence there should not be 
much difference in SMBH mass for galaxies with n >= 5. 
Above n ~ 12, increasing n to infinity has almost no effect 
on the profile shape, and hence the SMBH masses should 
all be the same. Some kind of asymptotic-like Mbh-n- func- 
tion is therefore in some sense demanded by the form of the 
Sersic R 1 ' 71 model. At the high-n end of the distribution, 
we do not believe that galaxies with n > 11.9 have smaller 
SMBH masses than those with n < 11.9, and we note that 
the log-quadratic Mbh-n relation (held fixed for n > 11.9) 
appears more logical than say a rising linear Mbh-n rela- 
tionship. In passing, we note that the highest SMBH mass 
which has been directly measured in a quasar using rever- 
beration mapping is only 2.6xl0 9 M Q (S50836+71, Kaspi 
et al. 2007) and is thus consistent with our predicted up- 
per 1-a range. The second highest quasar SMBH mass is 
O.9xlO 9 M (3C273). 

In this paper we employ the Mbh-n- relation in Eq[T]to 
derive the SMBH mass function using data from the Mil- 
lennium Galaxy Catalogue (MGC), which is described in 
Section 2. Preliminary results have been presented in Driver 
et al. (2006b, 2006c). Using the bulge-disc decompositions of 
the brightest 10 4 MGC galaxies (Allen et al. 2006), in Sec- 
tion 3 we construct the SMBH mass function. In Section 4 
we compare our results with previous efforts to measure the 
SMBH mass density, /Obh.o, using other means. A summary 
of our analysis is provided in Section 5. 

Throughout this paper, unless specified otherwise, we 
use Ov = 0.7, Q M = 0.3 and h 70 = H /(70 km s _1 Mpc" 1 ). 
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2 THE MILLENNIUM GALAXY CATALOGUE, 
AND OUR SPHEROID SAMPLE 

The MGC is a medium-deep (-Bmgc = 24 magjfl imag- 
ing survey of the nearby universe with a median seeing 
of 1.27" and it has 96.1 per cent complete (99.8 per cent 
for Bmgc < 19.2 mag) redshift information for the MGC- 
B RIGHT sample of 10,095 objects with Bmgc < 20 mag 
(Driver et al. 2005). The imaging data was acquired with the 
2.5 m Isaac Newton Telescope which surveyed 37.5 square 
degree^] in a 35 arcmin wide strip along the equatorial sky 
from Wh to 14ft 50' (Liske et al. 2003). Each field was ob- 
served for 750 seconds through a Kitt Peak National Obser- 
vatory -B-band filter (4407 A). The survey reaches a depth 
of /xiimit = 26 mag arcsec~ 2 , with objects catalogued down 
to Bmgc = 24 mag. For comparison, the Sloan Digital Sky 
Survey Data Release 5 (Adelman-McCarthy et al. 2007) has 
a median (r-band) seeing of 1.4" and an effective exposure 
time of 54.1 seconds per band, leading to a p-band (4686A) 
magnitude limit of 22.2 mag (roughly on an AB system, and 
with B — g roughly one-third of a mag, Blanton & Roweis 
2007). 

The MGC redshift information has come from a number 
of sources, as previously detailed in Driver et al. (2005, their 
Table 1). The median redshift is 0.12 and sample selection 
effects are well understood (Driver et al. 2005; Liske et al. 
2006). 

Allen et al. (2006) have performed an i? 1/ '"-bulge plus 
exponential-disc decomposition, using GIM2D (Marleau & 
Simard 1998; Simard et al. 2002), for all 10,095 objects. In 
addition to the best-fitting bulge Sersic index n, GIM2D 
derives the (not necessarily symmetric) upper and lower un- 
certainty, Sn, on the Sersic index. In FigO we show these 
uncertainties in n as a function of n. The distribution is 
such that 68 per cent of the galaxies have an error on n 
of less than ~20 per cent. Repeat observations, under dif- 
ferent seeing conditions and on different chips of the wide- 
field camera, exist for 682 galaxies. Fig. 15 in Allen et al. 
(2006) shows the ability of GIM2D to consistently recover 
the Sersic index of the spheroid component. The mean offset 
and standard deviation in the quantity A log n from repeat 
observations was reported there to be —0.002 and 0.132 dex, 
respectively. However, that distribution has longer tails than 
expected from a Gaussian. This is due to 11 per cent of the 
objects whose galaxy 'type' (see Allen et al. 2006) was in dis- 
agreement. Excluding these object^], the half-width of the 
68 percentile is 0.0729 dex, which corresponds to an 18 per 
cent mismatch in the value of n. This figure agrees well with 
the formal GIM2D error observed in Fig[3] and also with the 
20 per cent uncertainty used in Graham & Driver (2007a) 
and commonly reported in the literature (e.g. MacArthur, 

5 -Bmgc i s the Galactic extinction corrected, SExtractor (BEST, 
Vega) apparent magnitude (Liske et al. 2003). 

6 The actual usable area of sky reduced to 30.88 square degrees 
after excluding the "bad" regions (around bright stars, diffraction 
spikes, CCD defects, CCD gaps, CCD edges, vignetted corner, 
etc.). 

7 The half- width of the 68 percentile of the full distribution (i.e., 
all 682 galaxies) is 0.108 dex, which translates to a 28 per cent 
mismatch in the value of n, slightly higher than the formal GIM2D 
(lcr) error. This means that ~11 per cent of the data may have a 
larger uncertainty than is assigned by GIM2D. 



Courteau & Holtzman 2003). In this paper we adopted the 
GIM2D-derived values for both n and 5n, as given in the cat- 
alogue 'mgc_gim2d' which is publicly available at the MGC 
websitfl 

From the 10,095 galaxies in the MGC with -Bmgc < 20 
mag, we restrict the sample to those 7,745 objects with 
0.013 < z < 0.18. Many of these are disc-only systems 
and therefore rejected as potential black hole hosts. In pass- 
ing we note that, once the attenuating effects of dust have 
been dealt with, the MGC displays a uniformly flat distri- 
bution in cos(i), where i is the inclination of the disk, as one 
would expect for a uniformly distributed sample of galax- 
ies (Driver et al. 2007b, their figure 5). We further refine 
our galaxy sample by imposing the requirement that both 
the bulges and the elliptical galaxies, collectively referred to 
as 'spheroids', have half-light radii greater than 0.333 arc- 
sec (1 pixel) and that the bulge-to-total luminosity ratio 
(B/T) is greater than 0.01 (based on the lower values ob- 
served by Graham 2001, his figure 15). This helped avoid 
bright nuclear components such as star clusters that may 
have been fitted with the R 1/n model in GIM2D. We also 
required that the absolute spheroid luminosity be brighter 
than —18 B-mag (discussed further in Section |3.2.4|I . Fi- 
nally, we required that the galaxies core colour be red, such 
that (u — r) corc > 2.00 mag, denoting the transition in the 
colour bimodality for the MGC (Driver et al. 2006a) . We re- 
fer to the red sample as 'Sample 3'. However in Section \5. 2 1 
we show that the effect of including galaxies with blue cores 
(Sample 1 and 2) does not significantly alter our results on 
the SMBH mass density. 

We additionally construct two mutually exclusive sub- 
samples, which we label 'early-types' and 'late-types'. The 
distinction is based on whether a galaxy's B/T ratio is 
greater than or less than 0.4. Our choice of 0.4 is lower 
than values of 0.5 or 0.6 which have often been used in the 
past. Figure [4] shows the B/T ratio for ~3k MGC galaxies 
brighter than -Bmgc = 19 mag and which satisfy the above 
criteria and which have also been classified by eye into three 
morphological bins (Driver et al. 2006a): early- type galaxies 
(E/S0), early- type spiral galaxies (Sabc) and late- type spi- 
ral galaxies (Sd/Irr). If one was to use a cut at B/T = 0.5 or 
0.6 to identify the early-type galaxies, then one would miss 
roughly half of them. 



2.1 Colour correction 

Before proceeding, we note that Eq. [T] was constructed in 
the -R-band, while the MGC galaxies have been imaged in 
the B-band. The presence of radial, B — R colour gradients 
(e.g., La Barbera et al. 2005, and references therein) may 
therefore result in different values for the Sersic index n in 
the two bands. 

In general, colour gradients are known to be fairly small 
in observations of local, early-type galaxies more luminous 
than ~ — 17 -B-mag (e.g., Peletier et al. 1990; Taylor et al. 
2005, and references therein). The SDSS-VAGC (Blanton et 
al. 2005) has Sersic indices in the ugriz passbands in their 
low-redshift catalogue. Although these indices are somewhat 
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different to ours, in that they are derived from a single R ' n - 
galaxy model rather than an Bl^-bulge plus exponential- 
disc model, we should be able to get some insight from these 
data for the early-type galaxies, or at least the (disc-less) 
elliptical galaxies. An analysis of Fig[5] which shows the dif- 
ference between the SDSS r- and g-band Sersic index, re- 
veals that the median value of log(n r /n 9 ) for galaxies with 
n r > 2.0 — i.e. predominantly the early-type galaxies — is 
only 0.003 dex (after removal of the few obvious outliers with 
absolute values greater than 0.3 dex). We therefore apply 
no correction to our B-band Sersic indices of the early-type 
(B/T > 0.4) galaxies. 

The single B 1//n -galaxy models that have been applied 
to the SDSS data are not suitable for quantifying possible 
changes, with wavelength, to the Sersic indices of bulges in 
late-type galaxies. Instead, we use the B- and B-band Sersic 
indices from the bulges of 86 disc galaxies given by Graham 
(2003), 79 of which have indices in both bands. The average 
(± std.dev.) of the 79 values of \og{n R /n B ) is 0.09 (±0.15). 
Splitting the sample into Sa-Sb and Sc-Sd-Sm galaxies gave 
the same small offset of 0.09 dex for each grouping. This 
suggests that the difference in n is not a function of spi- 
ral galaxy type nor bulge size. In addition to radial stellar 
population gradients across the bulges, a plausible contrib- 
utor to this difference is dust. If the dust in spiral galaxies 
is more abundant at their centres, it will redden their cen- 
tres more than their outskirts, reducing the central parts of 
the B-band light profile relative to the B-band light profile 
and thereby yielding smaller Sersic indices for the bulge in 
the B-band. The similar B- and B-band Sersic indices for 
the early-type galaxies suggests that dust is not a signifi- 
cant issue in these systems. MacArthur et al. (2003) also 
provide B- and B-band Sersic indices for an independent 
sample of 47 and 43 late-type spiral galaxies, respectively 
— with 42 in common. In those instances where MacArthur 
et al. (2003) provided multiple Sersic indices for the same 
galaxy in the same passband, we averaged the logarithm of 
the Sersic indices. The average (± std.dev.) of the 42 val- 
ues of log(nfl/n_B) is 0.08 (±0.17), in good agreement with 
the data from Graham (2003). In deriving the SMBH mass 
function using Eq. [1] (established in the B-band), we will 
therefore apply a positive correction of 0.09 dex to the log- 
arithm of the B-band Sersic indices from the bulges of our 
MGC late-type galaxies. 

The presence of two distinct populations in FigO rather 
than one continuous distribution, suggests that the step-like 
correction we will apply to the Sersic indices of our early- 
and late-type galaxies (0.0 dex and 0.09 dex respectively) 
may be more appropriate than a continuous correction based 
on some parameter such as the B/T ratio. 



3 SMBH MASS FUNCTION AND SPACE 
DENSITY 

To derive the SMBH mass function, we first modify the 
Sersic indices of the late- type systems by 0.09 dex to convert 
from the B-band to the B-band, in accord with the previ- 
ous Section. We then derive individual black hole masses for 
each spheroid using equation [T] For each black hole we de- 
termine an associated space-density weighting based on the 
MGC blue and red spheroid luminosity functions as derived 



Table 1. Supcrmassive black hole mass function data (corrected 
for Malmquist bias) for the full, early- and late-type galaxy sam- 
ple (Sample 3, see Section l5,2H shown in Fig. [5] The uncertainties 
given are the upper and lower quartiles (i.e. ±25 per cent) from 
extensive Monte Carlo realisation of the combined errors. 
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Note: number densities are scaled to per unit logM^h interval and 
not per 0.251ogMt,h interval. 

in Driver et al. (2007a). The weight is the space density, 
4>(L), of the appropriate spheroid type (red or blue, divided 
at (it — r) C orc = 2.0 mag) in the specified luminosity interval, 
L, divided by the number of galaxies which contributed to 
that interval, N(L). The red and blue spheroid luminosity 
functions of Driver et al. (2007a) were derived via 'vanilla' 
step-wise maximum likelihood (see Efstathiou, Ellis & Pe- 
terson 1988) with K-corrections and e-corrections as defined 
in Driver et al. (2007a). The red and blue spheroid functions 
show markedly different forms and the justification for the 
segregation into two colour types is also given in Driver et 
al. (2007a). 

The SMBH mass functions are then derived by sum- 
ming the distribution of black hole mass times weights, i.e., 
<XM bh ) = £W(L)M bh , where W{L) = 4>{L)/N{L). This 
was constructed for black holes derived from all galaxies, 
early-types only (B/T > 0.4) and late-types only (B/T < 
0.4). See Fig[6]for a graphic representations and Table [1] for 
a tabulated version of these distributions. 

There are a number of sources of potential error. We 
model these via Monte Carlo simulations, following both in- 
dividually and collectively the uncertainties on: the param- 
eters defining equation [T] the Sersic indices; the luminosity 
function weights; and the Malmquist biafl Errors not mod- 
elled at this stage because they are considered secondary 
are uncertainties in the: K-corrections; e-corrections; dust 
attenuation; spectroscopic incompleteness; photometric er- 
rors effecting the magnitudes and hence weights; luminosity 
dependent cosmic variance; and the choice of colour cut in 

9 Here we define Malmquist bias to be the systematic bias in our 
final measurements due to the presence of errors in our data. 
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the derived red and blue spheroid luminosity functions. Note 
that the global cosmic variance for the MGC survey area was 
derived in Driver et al. (2005) by comparison to mock cat- 
alogues produced by the Durham grourj"! this amounts to 
an overall uncertainty of 6 per cent in all 4 x 10, 001 Monte 
Carlo simulations. 

The Monte Carlo simulations consist of repeating the 
above analysis 10,001 times with the input values per- 
turbed by the listed errors individually and collectively, i.e., 
4 x 10, 001 Monte Carlo simulations in all. All error dis- 
tributions are assumed to be Gaussian. The impact of the 
error(s) is(are) then assessed by comparing the median and 
standard deviations of the derived distributions to the orig- 
inal estimate. 

Table [2] summarises the cumulative sum of the SMBH 
mass functions, Pbh.o, with error estimates from each of the 
above sources. In all cases we see that the dominant error 
is from the uncertainty in equation [1] This is perhaps not 
surprising since equation Q] is based on a quadratic fit to 
only 27 systems for which both credible black hole masses 
and Sersic indices exist. Improvements in this method will 
therefore come from a larger calibration sample (similarly 
for the M — a and M — L relations) and this should be 
attainable with next generation facilities. Perhaps surprising 
is that the errors in the Sersic distribution have relatively 
little impact. In part this is indicative of the size of the 
sample but also helped by the quadratic nature of the M — n 
relation preventing high n values leading to unreasonably 
large Mbh values. The errors from the luminosity function 
(i.e., statistics and cosmic variance) appear negligible, with 
cosmic variance the most significant. 

Finally, we correct for the Malmquist bias in a rather 
straightforward manner. This is achieved by measuring pbh.o 
assuming no errors and then remeasuring it but allowing 
all of the above errors to be perturbed simultaneously via 
Monte Carlo simulations. The difference between the two 
estimated values (raw and the median of the Monte Carlo 
simulations) provides a crude measurement of the system- 
atic offset caused by the errors inherent in our data. To 
provide Malmquist bias corrected values we then subtract 
this systematic from our original values (~ 1 to 3 per cent 
downward correction). These final Malmquist bias corrected 
values are those shown in Tables [T] and [5] 

FigHJ] shows the resulting black hole mass functions, 
where the error bars indicate the combined errors derived 
from the Monte Carlo analysis but excluding the cosmic 
variance error (as this is a hidden systematic and not a 
random error). Our resulting SMBH mass functions (and 
number density and mass density) depend on the Hubble 
constant. As was shown in Graham & Driver (2007a), the 
Mbh-n relation is independent of the Hubble constant. This 
is because the Sersic index n does not depend on galaxy dis- 
tance, and while the SMBH mass does depend on distance, 
the overwhelming majority of galaxies that were used to 
construct the Mbh~n relation had their distances obtained 
by Tonry et al. (2001) using surface brightness fluctuations 
- a technique that provides distances without assuming 
some Hubble constant. The mass function dependence on 
the Hubble constant arises from the 1/Vmax weighting given 
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Table 2. The space density of matter in supermassivc black holes. 
The errors are 1-cr values, after excluding 3-cr outliers. The densi- 
ties given in Table[T] rather than the fitted empirical models, have 
been integrated down to SMBH masses of 10 6 M Q . The final col- 
umn shows the density normalised against the critical density. The 
reason why the density varies with h? and not h 2 is explained in 
Section 3. Factoring in the intrinsic scatter, A, from the M^^-n 
log-quadratic relation, the numbers in this Table should be in- 
creased by exp[(AlnlO) 2 /2], which equals 1.09 if A = 0.18 dcx 
(Graham & Driver 2007a). 







No. 


Pbh.ql ±<5^ ±5£ ±5^ ±8^_ 


^bh 






Bulges 


ft| 10 5 M Q Mpc- 3 








Early 


- and Latc-typc (B/T > 0.01) 




Sample 


1 


1769 


4.87 ± 1.84 ± 0.07 ± 0.01 ± 0.29 


3.6±1.4 


Sample 


2 


1676 


4.53 ± 1.63 ± 0.07 ± 0.01 ± 0.27 


3.3±1.2 


Sample 


3 


1543 


4.41 ± 1.65 ± 0.07 ± 0.01 ± 0.26 


3.2±1.2 








Early-type (B/T > 0.4) 




Sample 


1 


1539 


3.78 ± 1.28 ± 0.06 ± 0.01 ± 0.23 


2.8±1.0 


Sample 


2 


1485 


3.57 ± 1.21 ± 0.06 ± 0.01 ± 0.21 


2.6±0.9 


Sample 


3 


1352 


3.46 ± 1.14 ± 0.06 ± 0.01 ± 0.21 


2.5±0.9 






Late- type (0.01 < B/T < 0.4) 




Sample 


1 


230 


1.14 ± 0.56 ± 0.04 ± 0.01 ± 0.07 


0.8±0.4 


Sample 


2 


191 


0.95 ± 0.49 ± 0.03 ± 0.01 ± 0.06 


0.7±0.4 


Sample 


3 


191 


as above 





a Malmquist bias corrected. 

^ Monte Carlo simulation of the errors in the M^-xi relation 
assuming a Gaussian error distribution. 

c Monte Carlo simulation of the error in n as specified by GIM2D 
and assuming a Gaussian error distribution. 

^ Monte Carlo simulations of the error in the individual galaxy 
weights assuming a Gaussian error distribution. 
e MGC global cosmic variance of 6 per cent for the effective ~ 30 
sq degree region with 0.013 < z < 0.18. 

to each galaxy; it is this volume term which introduces an 
h 3 dependence. 



3.1 Parameterisation of the SMBH mass function 

We have fitted two empirical models to our SMBH mass 
functions over the mass range 10 6 < M bh /M© < 10 9 . The 
first is a mild variation of the commonly-used 3-parameter 
Schechter (1976) function, and is given by 



0(M bh ) = 0, 



exp 



_ / M bh 
V M» 



(2) 



where (j>(M*) = <f) t (per unit dlog(Af b h) per Mpc 3 ). The 
turnover of the mass function and the maximum density 
occur at the SMBH mass 



M max = (a+ 1)M», a + l>0 
where the associated maximum density is 
0max = <A»(a + l) a+1 exp [-a] . 



(3) 



(4) 



The logarithmic slope at the low-mass end of the mass 
function is given by the exponent 1 + a; a value of a = 
— 1 therefore corresponds to a flat distribution, and larger 
values correspond to a decreasing function as the SMBH 
mass decreases. For the early-type galaxy samples the slope 
(1+a) is ~1 (see Table[3]). For the (early+late)-type samples 
the slope is approximately two-thirds; the shallower decline 
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is due to the contribution of SMBHs in late-type galaxies. 
These, and the other, best-fitting parameters are reported 
in Table [3] and the fits themselves are shown in Fig[6] 

We explored the suitability of a second model having 
an additional fourth parameter such that the mass term in 
the exponential of equation [2] is raised to the power of j3 to 
give 



0(M bh ) = 4>* 



Mbh 
Ah 



exp 



(Mbh\> 
\ M* J 



(5) 



(Aller & Richstone, their equation 10). Our data, how- 
ever, did not justify the need for this additional parameter. 
For the (early+late)-type samples, the value of j3 equalled 
1.0±0.1. For the early-type sample, it ranged from 0.3 to 
0.5, but did not give rise to significantly better fits. For this 
reason we do not show the fits or the parameters from this 
model. (Neither equations[2]nor[5]provided an acceptable de- 
scription to the SMBH mass function in our late-type galaxy 
sample.) 



3.2 Integrating the SMBH mass function 

8.2.1 SMBH number density 

The total SMBH number density can be obtained by inte- 
grating equation [5] with respect to log(Mbh) (because <j) is 
expressed in units of /i 3 Mpc" 3 per decade in SMBH mass). 
Integrating over 10 6 < Mbh/M© < 10 9 , the number density 
is given by the expression 



l°g(M bh /M )=9 



<j>(M hh ) d log Mbh 



A/ bh = 10 a M o 



</>{Mbh) 

M bh =io«M 8 ln(10)Af bh 
^e 1 r ( . , 10 9 M, 



dM b 



ln(10) 



7 [a + 1. 



M, 



7 [a + 1 



10 6 M ( 



Ah 



((».) 



where 7(0, x) is the incomplete gamma function (e.g., Press 
et al. 1992) defined by 



7(0,0:) = 



(7) 



For the early- and late- type galaxy values of <j>* , A'h and 
q (given in Table the number density of SMBHs with 
masses between one million and one billion solar masses is 
2.3 x 10~ 3 /i 3 Mpc -3 . For the early-type galaxies, the num- 
ber density is 2.1 x 10" 3 /i 3 Mpc" 3 . 

These values are an order of magnitude smaller than 
those reported in Shankar et al. (2004). The difference can 
be attributed to their rising (Mbh--^)-derived SMBH mass 
function (for all galaxy types) as one moves towards lower- 
masses. In contrast, our (Mbh-n)-derived mass function 
shows the opposite behaviour at the low-mass end. 



3.2.2 SMBH mass density 

In so far as equation [2] represents the mass function over the 
SMBH mass range 10 6 < M bh /M @ < 10 9 , the SMBH mass 



density for SMBHs having such masses can be obtained from 

/■log(M bh /M )=9 

pbh = / (j)(M bh )M bh dlogMbh 

•Aog(A/ bh /A/ )=6 



ln(10) 



7 I a + 2. 



1O 9 M 
Ah 



-7 [a + 2, 



1O 6 M 

A'h 



(8) 



Using the best-fitting parameters in Table [3] we obtain 
Pbh = 4.3 x 10 5 /i 3 M Q Mpc" 3 (all-types), and p hh = 
4.0 x 10 5 /i 3 o-Mq Mpc" 3 (early- types). Integrating the mass 
function over all masses increases the above values to 4.8 x 
10 5 /i 3 M© Mpc" 3 and 4.2 x 10 5 /i 3 M© Mpc" 3 , respectively. 

We have, however, opted not to use the above formula, 
but to instead acquire the SMBH mass densities directly 
from our data, rather than from the fitted model. 

Computing the local mass density, and integrating down 
to SMBH masses of 10 6 A4q, we obtain phh,eariy-type = (3.5± 
1.2) x 10 5 /i? Af Q Mpc" 3 while pbh,iatc-typo = (1.0 ± 0.5) x 
l0 5 h,7 M Q Mpc" 3 . These results are presented in Tabled 

The above densities correspond to f2bh,eariy-ty P e = 

Pbh,carly-typc/pcrit = (2.5 ± 0.9) X 10" 6 /l 7 0, and 

^bh,iato-t yP o = (0.7 ± 0.4) x 10" e /i 7 o (see Table 0, where 
Pcrit = 3Hq/8ttG is the critical density for flat spacj"!. 
For reference, using an independent technique, Fukugita & 
Peebles (2004, their equation 75) give 2. 5±f l x 10" 6 for the 
SMBHs in early-type galaxies, and 1.31q^ x 10" 6 for the 
SMBHs in late-type galaxies. 

In deriving the (linear) SMBH mass density, Pbh.o, (i- e -, 
not log Pbh.o) from the convolution of the distribution func- 
tion of n with the log Mbh-log n relation, one needs to allow 
for the intrinsic scatter in this log — log correlation. If there 
is a Gaussian distribution of intrinsic scatter in the SMBH 
mass at any fixed logn, with a standard deviation which is 
independent of logn and equal to A, then the SMBH mass 
density should be increased by the factor exp[(A In 10) 2 /2] 
(Yu & Tremaine 2002, their equation 12). Such a Gaus- 
sian distribution, however, is not true for the scatter of 
points about the log-quadratic M bb -n relation shown in 
Figtll the scatter is clearly greater at the high-mass end. 
In fact, removal of the two highest mass SMBHs results in a 
log-quadratic relation consistent with zero intrinsic scatter 
(Graham & Driver 2007a). It is thus questionable whether 
one should apply this multiplicative term and because of 
this uncertainty, in what follows, we have not. However, one 
should note that if there were in fact no measurement er- 
rors and A = 0.31 dex (the total absolute scatter about the 
log Afbh-log n relation in FigfTJ , then the corrective factor to 
apply to pbh,o would be 1.29. If A = 0.18 dex (the intrinsic 
scatter after factoring in the suspected measurement errors), 
then the multiplicative factor drops to 1.09. If there is no 
intrinsic scatter, then the multiplicative factor is simply 1. 



3. 2. 3 SMBH baryon fraction 

The above numbers can be expressed in terms of the baryon 
fraction of the universe. Already a picture is emerging 
in which dark accretion onto SMBHs may be negligible. 
Shankar et al. (2004) claim that the local SMBH mass func- 
tion can be accounted for from mass accreted by X-ray se- 



11 For H = 70 km s" 1 Mpc" 1 , p crit ~ 1.36 x 10 n M Q Mpc" 
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Table 3. Best-fitting parameters from the empirical SMBH mass function given in 
equation [2] The number density, </>„, is per decade in SMBH mass. The late-type 
galaxy sample is not shown here due to the poor match between the empirical model 
and the data (see Figj6j . Sample 2 excludes bulges in disc galaxies if the core colour is 
bluer than (u — r) c = 2.00 mag. Sample 3 (our primary sample) excludes blue bulges 
and blue elliptical galaxies if (u — r) c < 2.00 mag (see Section l5.2l l. 



Data sample 




log(M,/M ) 


a 




h'^ Q Mpc~ 3 dex -1 






Early- and Late-type, (B/T > 0.01): 








Sample 1: no colour cut 


-2.76 


8.45 


-0.32 


Sample 2: (u — r) c [bulges] > 2.00 mag 


-2.76 


8.43 


-0.29 


Sample 3: (u — r) c [all] > 2.00 mag 


-2.81 


8.46 


-0.30 


Early-type, (B/T > 0.4): 








Sample 1: no colour cut 


-2.68 


8.26 


0.07 


Sample 2: (u — r) c [bulges] > 2.00 mag 


-2.67 


8.25 


0.10 


Sample 3: (u — r) c [all] > 2.00 mag 


-2.74 


8.29 


0.00 



lected AGN, i.e. no significant 'dark' accretion is required 
for SMBH growth (see also Cao 2007). If correct, this would 
imply that the growth of SMBHs from the accretion of mas- 
sive black hole remnants of Population III stars may be small 
(Madau & Rees 2001; Islam, Taylor & Joseph 2003; Wyithe 
& Loeb 2004). Moreover, baryonic-fuelling is consistent with 
a picture that links the BH mass to the host spheroid's bary- 
onic, or at least stellar, properties, such as luminosity (e.g., 
Erwin, Graham & Caon 2002; McLure & Dunlop 2002; Mar- 
coni & Hunt 2003 and references therein), concentration 
(Graham et al. 2001; Graham & Driver 2007a) and mass 
(Magorrian et al. 1998; Haring & Rix 2004). Assuming the 
above standard picture in which SMBHs form via the ac- 
cretion of baryons (Blandford 2004), then in terms of the 
invariant baryon fraction of the total mass-energy density, 
such that fibaryon = 0.0453/iyQ 2 (Tegmark et al. 2006, their 
Table 2, row 3; see also Blake et al. 2007), SMBHs at the 
centres of galaxies today contain (0.007 ± 0.003)/i7 per cent 
of the universe 's baryon inventory. 

Letting M sp h C roid denote the stellar mass of a spheroid, 
Haring & Rix (2004) have shown that Mbh/A^sphcroid = 
0.0009 when Af bh = lO 6 Af and = 0.0019 when M bh = 
10 9 Mq, confirming the results given in Merritt & Ferrarese 
(2001, their Section 3; see also Laor 2001). Taking the av- 
erage (or using the peak in our SMBH mass function at 
log(Mbh) = 8.3, see Fig|6j one obtains Mbh equals 0.14 per 
cent of Af spheroid (or 0.16 per cent). Dividing fibh/^baryon 
by this value one obtains fi sp hcroid ~ 6 per cent (or ~5 per 
cent) of fibaryon ■ This is of course just a rough estimate, and 
a more precise value obtained using the actual spheroid lu- 
minosity function will be presented in Driver et al. (2007a, 
2007b). 



3.2.4 Comments on the low mass end 

Figure [7] shows how the SMBH space density is built up 
as one includes (intrinsically) fainter spheroids, for our 
three galaxy bins. As was noted in Section 2, we excluded 
spheroids fainter than Mb — — 18 mag in our analysis. For 
the early-type galaxy sample, Figure [7] reveals that integrat- 
ing down to an absolute magnitude of —16 B-mag gives a 
SMBH mass density that is consistent (at the 2-cr level) with 
the result obtained using only those spheroids brighter than 
— 18 B-mag. We do however caution that the SMBH masses 



are less reliable in this regime, and as such we do not wish 
to place too much emphasis on this result. 

It is also worth noting that the central massive object 
in spheroids fainter than Mb ~ —18 mag is often a nuclear 
star cluster (Ferrarese et al. 2006, left panel of their fig- 
ure 2; Balcells et al. 2007). It is therefore questionable as to 
whether or not one should include such faint spheroids in an 
analysis of this kind. Moreover, if they are not individually 
modelled, the presence of additional nuclear star clusters 
can bias the Sersic R 1 ^" fit to give spuriously high values of 
the index n. The effect is to over-estimate the contribution 
from SMBHs in faint spheroids. Given the high frequency 
of nucleated bulges in spiral galaxies (e.g., Carollo, Stiavelli, 
& Mack 1998; Balcells et al. 2003), there is thus reason to 
doubt the rising SMBH space density curve shown in Fig[7] 
for the late-type galaxies. 



4 COMPARISON WITH PAST ESTIMATES 

McLure & Dunlop (2004) provided two estimates of the 
SMBH mass function (in early-type galaxies) using two tech- 
niques. At the low-mass end (< 10 s Mq) their estimates do 
not agree; the mass function they obtained using the Mbh - 
L relation gave noticeably higher SMBH number densities 
than they obtained using the M^-a relation (Fig(8] lower 
panel). Both the method and data that we have used is dif- 
ferent to that used in McLure & Dunlop (2004) and thus 
provides an independent check on the SMBH mass func- 
tion. Our analysis, using the Mbh-^ relation, shows a mass 
function which declines with decreasing SMBH masses that 
are less than 10 8 Mq, and therefore better matches McLure 
& Dunlop's mass function constructed using the Mbh-c re- 
lation rather than the M^-L relation. 

In plotting the mass functions from McLure & Dunlop 
(2004) in FiglS] we have multiplied their (M b h-£)-derived 
masses by 1.70. This increase stems from the conversion of 
their B-band Mhh~L relation to the if-band. Starting from 
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equation 6 in McLure & Dunlop (20020 , one has 

log(M bh /M ) = -0.5M R - 2.91 
= -0.5[Mx + 2.7 - M K ,@ + 3.28 - 51og(70/50)] - 2.91 

= 1.25 log(L K /L K;@ )- 5.53 (9) 

Following McLure & Dunlop (2002), we have assumed an 
average R — K colour of 2.7. We have used Mk,q = 3.28 mag 
(Binney & Merrifield 1998). McLure & Dunlop (2002) used 
the SMBH masses derived using Tonry et al.'s (2001) surface 
brightness fluctuation distances — which are independent 
of the Hubble constant. However they did not use these h- 
independent distances in deriving the absolute magnitude of 
the bulges, but used redshift derived distances and Ho = 50 
km s -1 Mpc -1 . Equation[9]gives SMBH masses that are 0.23 
dex larger than equation 1 in McLure & Dunlop (2004), and 
hence the factor of 1.70 (see Graham & Driver 2007b for 
more details). 

In the lower panel of Fig[S] one can see that our analysis 
suggests there may be a greater (up to a factor of ~2) num- 
ber density of SMBHs with masses log(Mbh/MQ) ~ 8.5 than 
reported in McLure & Dunlop (2004). At smaller masses, 
8 > log(Mbh/MQ) > 6, our mass functions roughly display 
the same decline with decreasing SMBH mass as McLure & 
Dunlop's analysis based on the M^~a relation. Such a de- 
clining SMBH mass function (for early-type galaxies) is also 
in qualitative agreement with Shankar et al.'s (2004) mass 
function derived using a bivariate velocity dispersion distri- 
bution. In passing we note that we have increased Shankar 
et al.'s SMBH masses by 14 per cent, after adjusting for 
the correct dependence on the Hubble constant (see Gra- 
ham & Driver 2007b). On the other hand, our data dis- 
agree with mass functions that rise monotonically (which 
includes asymptotically) towards lower-mass SMBHs, such 
as the (Mbh-i)-derived mass functions in McLure & Dunlop 
(2002), Aller & Richstone (2002) and Shankar et al. (2004). 

The upper panel of Fig(5] reveals a clear disagreement, 
at masses below ~ 10 s Mq, between our (total) mass func- 
tion and that derived from Mbh~L relations. We are not 
aware of any (Mbh _ o")-derived mass function for late-type 
galaxies, and therefore we do not show any such curve in 
this panel. Some of the mismatch from the (M b h-I/)-derived 
SMBH masses may stem from the methods used to assign 
bulge flux in lenticular and spiral galaxies (see Graham & 
Driver 2007b and Graham 2007). 

Over the past few years a number of papers have con- 
structed the local SMBH mass function and estimated the 
local SMBH mass density, pbh.o- We have listed several of 
these in Table [4] It can be seen that our estimates of pbh,o 
agree particularly well with the estimate in Fukugita & Pee- 
bles (2004), and is consistent (within the 1-a error bounds) 
with several other recent studies. It is worth noting, however, 
that our estimates are a factor of ~2 greater than reported 
in Wyithe (2006) and Aller & Richstone (2002). Due to the 
declining number density with decreasing SMBH mass, our 

12 Transformation of equation 5 in McLure & Dunlop (2002) 
would be considerably more difficult because it was constructed 
from an inactive plus active galaxy sample, a fraction of which 
have distances that depend on the Hubble constant and also on 
the adopted cosmology given that some of the active galaxies have 
redshifts which extend out to ~ 0.5. 



adopted integration across the complete mass spectrum of 
SMBHs, rather than truncating at 1O 6 M0, does not signifi- 
cantly increase our estimate of Pbh.o- 

In an effort to try and understand why our SMBH mass 
densities may be slightly higher than reported by some, 
Graham & Driver (2007b) has examined two representative 
studies; one which used the M^-a relation (Aller & Rich- 
stone 2002) and the other the Mbh~L relation (Shankar et al. 
2004). Graham & Driver identified a number of corrections 
which can be made and showed that the mass density from 
Aller & Richstone (2002) is more than a factor of ~ 2 too 
low. As with the analysis by McLure & Dunlop, which was 
a factor of 1.7 too low, the reason is because of over-looked 
dependencies on the Hubble constant. These have been cor- 
rected in Figure [SJ but not incorporated into Table [3] 



5 POTENTIAL BIASES IN OUR DATA 
5.1 A Linear Mbh~ n relation 

While we believe the quadratic relation between Mbh and n 
(equation]!} is the optimal expression to use when predicting 
Mbh from n, we have explored how pbh.o would change if we 
adopted the linear relation presented in Graham & Driver 
(2007a, their equation 2). Although a linear fit to the data 
in Figure [1] results in predictions of larger SMBH masses at 
the high-n end, over the range ~2 < n < ~8 it actually 
predicts lower SMBH masses (see Figure 3 in Graham & 
Driver 2007a). Because the bulk of the MGC data falls in 
this interval (see Figure [o r )P^I. use of the linear M b h-n rela- 
tion results in a 34 per cent lower SMBH mass density. This 
value is consistent, i.e. within the 1-a uncertainty, with our 
estimate p b h,o = (4.4 ± 1.7) x 10 5 fr 3 M Q Mpc -3 . However, 
as noted in Graham & Driver (2007a), this linear relation 
has an intrinsic scatter of 0.31 dex, and so the estimate of 
Pbh.o should be increased by 29 per cent, to give a value of 
3.8 x 10 5 /i 3 o M Mpc" 3 . 

We have additionally removed the three high-n data 
points from Figure[T]and obtained a new linear relation with 
the remaining 24 data points. Doing so results in a relation 
with zero intrinsic scatter and a mass density that is 43 
per cent higher than obtained with the quadratic relation. 
Specifically, we obtain a value of (6.3 ± 2.5) x 10 5 /i? M s 
Mpc -3 , which is only marginally higher than our upper 1-a 
value of 6.1 x 10 5 /i 3 oM Q Mpc -3 . We also find Pbh.early-type = 
(4.0 ± 1.2) x 10 5 h%M Q Mpc -3 and p b h,iatc-ty pc = (2.3 ± 
1.3) x W 5 h% M Q Mpc -3 . 

The associated SMBH mass functions are shown in Fig- 
ure[5]and given in Tabled We caution that we do not believe 

13 We speculate that the apparent shortage of Sersic indices 
greater than ~8 in the MGC data (compared to Figure [TJ may 
arise from GIM2D's difficulty in deriving such large values while 
producing a galaxy magnitude that is close to the value obtained 
using SExtractor (due to the large wings of high-n profiles). Con- 
firmation of this speculation would require simulations beyond 
the scope of this paper. We therefore offer it as nothing more 
than speculation, but do note that if correct, a variety of plausi- 
ble stretches to the MGC values of n that we tested (to recover the 
possible true distribution) resulted in a ~20 per cent increase to 
the value of pbh o (i- e -i smaller than our current 1-a uncertainties 
on this value). 
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Table 4. Local SMBH mass density estimates. The difference between Samples 1 and 3 (see 
Section IBT2t in our study is that the latter excludes galaxies with (u — r) COTC < 2.00 mag, such 
as blue pseudo-bulges. The factor h 3 = [Ho/(70 km s _1 Mpc -1 )] 3 is appropriate for our study 
because the M^-n relation is independent of the Hubble constant. The majority of the densities 
from other papers have been transformed to Hq = 70 km s" 1 Mpc -1 using b? rather than /i 3 , 
as indicated in each paper. However, as discussed in Graham & Driver (2007b), this may not 
always be appropriate. 



Study 


p hh . (E/S0) 


Pbh,o (Sp) 


Pbh.o (total) 




h^ W 5 M Q Mpc" 3 


h^ W 5 M Q Mpc" 3 


h% 10 5 M Q Mpc" 3 


This study (Sample 1) 


(3.8±1.3)/i 7 o 


(1.1 ± 0.6)h 7Q 


(4.9 ± l.9)hro 


This study (Sample 3) 


(3.5 ± 1.2)h 70 


(1.0±0.5)/i 70 


(4.4 ± 1.7)h m 


Wyithe (2006) 






2.28 ± 0.44 


Fukugita & Peebles (2004^ 




fl 7 +1 - 7 )h- 1 


(5.ll?:2)hro 


Marconi et al. (2004) 


3.3 


1.3 




Shankar et al. (2004^ 


o 1+0.9 

^•-'■-o.s 


— 0.5 




McLure & Dunlop (2004) 


2.8 ±0.4 






Wyithe & Loeb (2003) 






9 9 +3.9 


Aller & Richstone (2002^ 


1.8 ±0.6 


0.6 ±0.5 


2.4 ±0.8 


Yu & Tremaine (2002)£ 


2.0 ±0.2 


0.9 ±0.2 


2.9 ±0.4 


Merritt & Ferrarese (2001)^ 






4.6/ifo 1 


Salucci et al. (1999) 


6.2 


2.0 


8.2 



a See their equation 75. 

^ Based on their (Mbh-i)-derived mass function, and in agreement with their (Mt,h-0-)-derived 
values. 

c Taken from their Table 2. 

^ Based on their (M^— <r)-derived mass function. 
e See also Ferrarese (2002). 



these are more accurate than those provided in Table [T] but 
that they provide an (extreme) upper limit to the SMBH 
mass function and space density. As was noted in Section 2, 
the evidence for a curved M^-n relation is inherent in the 
low-n data. Excluding the five highest n galaxies results in 
a quadratic relation that is fully consistent with equation [1] 
Removing the two galaxies with the highest SMBH masses 
(one of which is suspected to be biased high: NGC 4486) 
produces the same result. The linear relation obtained after 
removing the three high n galaxies is thus entirely driven by 
two galaxies (NGC 4486 and NGC 4649). 



5.2 Blue Spheroids and Pseudobulges 

We divided each of the early-type (E/S0), late-type (Sp) 
and (early + late)-type galaxy bins into three (colour) sam- 
ples. Sample 1 has had no modifications to it, while Sam- 
ple 2 excludes bulges (in disc galaxies) if the (u — r) core 
colour is bluer than 2.00 mag. The core-colour is derived 
from Sloan Digital Sky Survey point-spread-function mag- 
nitudes which have been shown to maximise the colour bi- 
modality, see Driver et al. (2006a). This colour-cut may help 
avoid possible pseudo-bulges (e.g., Erwin et al. 2003; Kor- 
mendy & Kennicutt 2004; Kormendy & Fisher 2005) which 
may or may not have SMBHs. Sample 3 is further reduced 
to exclude all systems (i.e., bulges and elliptical galaxies) 
with (u ~ r) c < 2.00 mag (see Fig fTU|l . 

It is not clear whether supermassive black holes ex- 
ist in (every) pseudo-bulge. The tight correlation between 
the properties of dynamically hot spheroids and their black 
hole mass has been tied to their joint formation process. As 
pseudo-bulges are believed to have formed via a distinctly 



Table 5. Supermassive black hole mass function data (corrected 
for Malmquist bias) for the full, early- and late-type galaxy sam- 
ple (Sample 3, see Section 15.21 1 shown in Figure [9] The linear 
Mbh-n relation (Section [STTJ, obtained after excluding the three 
galaxies with the highest Sersic index, has been used here to pro- 
vide an upper estimate to the SMBH mass function. The uncer- 
tainties given are the upper and lower quartiles (i.e. ±25 per cent) 
from extensive Monte Carlo realisation of the combined errors. 
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Note: number densities are scaled to per unit logMbh interval and 
not per 0.251ogMj,h interval. 
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different formation process (i.e., from a redistribution of 
stars due to perpendicular-to-the-disc resonances with a bar) 
it is plausible that pseudo-bulges do not possess SMBHs. De- 
spite this conjecture, we note that bars can fuel AGN growth 
(e.g., Ann & Thakur 2005; Combes 2006; Ohta et al. 2006). 
Moreover, the Milky Way has a SMBH mass that appears 
consistent with both the Mbh~n and Mbh-o" relations and 
probably has a pseudo-bulge — suggested from its rotation 
(e.g., Minniti et al. 1992), flattening (e.g., Sharpies, Walker 
& Cropper 1990; Lopez-Corredoira et al. 2005), presence 
of a bar (e.g., Lopez-Corredoira et al. 2006) and peanut- 
shaped structure (e.g., Alard 2001, see also Debattista et al. 
2006), but see Zoccali et al. (2006) and Matteucci (2006). 
The situation is therefore ambiguous, which has resulted in 
our decision to create three subsamples. 

To unequivocally distinguish dynamically hot bulges 
from rotating pseudo-bulges requires kinematical informa- 
tion. In the absence of this, an effective method is to use 
colour, with the pseudo-bulges typically exhibiting blue 
colours comparable to a young (metal rich) disc population 
and the genuine bulges exhibiting red colours consistent with 
an old population. We do however note that in some cases 
the bulges are poorly resolved and the central colour is thus 
a blend of the bulge and disc, which might result in some 
real bulges erroneously falling blueward of our colour cut. 
However the fact that the colour bimodality is clearly evi- 
dent in our data suggests this is probably not a major issue, 
and we note that a positive correlation between bulge and 
disc colour is known to exist (e.g., Peletier & Balcells 1996). 

Including blue spheroids with (u — r) c < 2.00 mag, Sam- 
ple 1 gives densities of Pbh,cariy-t yP c = (3.8±1.3) x 10 5 /i 3 Mq 
Mpc -3 and p b h,iato-t ypc = (1.1 ± 0.6) x W 5 h% M Q Mpc" 3 
(Table I). 



6 SUMMARY 

We have used the log-quadratic relation between supermas- 
sive black hole mass and the host spheroid's Sersic index 
n (i.e. the M b h-fi relation, Graham & Driver 2007a) to 
predict the SMBH masses in galaxies from the Millennium 
Galaxy Catalogue (MGC) for which a two-dimensional, 
seeing-corrected, R 1 /™ -bulge plus exponential-disc decom- 
position has been performed (Allen et al. 2006). 

Applying the appropriate volume corrections, the num- 
ber density of SMBHs with masses greater than 10 6 Mq, 
and in spheroids more luminous than Mb = — 18 mag, is 
(2.3 ±0.1) x 10~ 3 /i? o Mpc -3 . 

We have constructed the local (0.013 < z < 0.18,2 = 
0.12) SMBH mass function using direct estimates of the 
Mbh predictor-quantity n (Fig|5J. Using every spheroid, 
the SMBH mass function appears well represented by a 3- 
parameter Schechter function having a logarithmic slope at 
the low-mass end of around two-thirds. Excluding the late- 
type galaxies, identified as systems with bulge-to-total ratios 
smaller than 0.4, the SMBH mass function was again well 
fit with a Schechter function, this time having a low-mass 
slope of 1, and a turnover mass at ~ 2 x 10 8 Mq. 

Integrating down to a SMBH mass of 10 6 Mq, the local 
mass density of SMBHs from early- and late-type galaxies 
combined is p b h,o = (4.4 ± 1.7) x 1O 5 /i 3 o M Mpc" 3 . This 
value is slightly greater than, but still consistent with, most 



values obtained via independent means over recent years 
(Table [4] but see Graham & Driver 2007b). Our mass den- 
sity estimate can be split into Pbh,eariy-t ypc = (3.5 ± 1.2) x 
W 5 h% M e Mpc" 3 and p b h,iatc-t yP c = (1.0±0.5) x 10 5 h% M Q 
Mpc -3 , based upon a B/T flux ratio cut at 0.4. 

In terms of the critical density of the universe, we obtain 

^bh, total = (3.2 ± 1.2) X 10 _6 /l70, Or fibh.carly-tyfjc = (2-5 ± 

0.9) X 10~%o and fibh,late-type = (0-7 ± 0.4) x W- 6 h 70 . 
Given that 4.5/if 2 per cent of the critical density is in the 
form of baryons (Tegmark et al. 2006), the above figures 
imply (0.007±0.003)/i 3 per cent of the total baryon content 
of the universe is locked up in SMBHs at the centres of 
galaxies. 
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Figure 1. Relationship between SMBH mass and host spheroid 
Sersic index. The data points were taken from Graham & Driver 
(2007a, their table 1) and the curve shows the best quadratic fit 
(equation 1). 
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Figure 4. Distribution of GIM2D-derived, bulge-to-total flux ra- 
tios (B/T) for those MGC galaxies with an (eyeball) morpho- 
logical classification (Driver et al. 2006a). Galaxies with a core 
colour bluer than (u — r) c = 2.00 mag have been excluded from 
the top panel. Spheroids fainter than Mg = —17 mag and effec- 
tive half-light radii smaller than 0.333 arcseconds (1 pixel) have 
been excluded from all panels. The dashed line at B/T = 0.4 de- 
notes our adopted divide between early- and late-type galaxies. 
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Figure 3. GIM2D-derived uncertainty in n, namely Sn (which 
need not be symmetrical about n), is converted to a fractional un- 
certainty and shown as a function of n for the full MGC spheroid 
sample. Sixty eight per cent of the data has an error on n of less 
than 20 per cent. The bottom-left envelope of points in the left 
panel trace the lower limit on the value of n (equal to 0.2) that 
we set in GIM2D (see Allen et al. 2006). A second envelope arises 
from the 0.09 dex colour-correction which is applied to the bulges 
of the late-type galaxies (see Section 12. U . 
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Figure 2. Left panel: Variant of Fig[T] after removing the five data points with the highest SMBH masses. The curved line corresponds 
to logM bll = a + log(n/3) + c[log(n/3)] 2 . Projections of the A\ 2 = 1.00,4.00 and 6.63 contours, shown in the central and right panels, 
onto each axis gives the 68.3, 95.4, and 99 per cent confidence interval on each of the parameters a, b and c. The optimal log-quadratic 
relation (shown in the left panel) is consistent with zero intrinsic scatter and reveals that the coefficient, c, in front of the quadratic term 
is inconsistent with a value of zero at the ~ 2.5 a level. Subsequently, one can conclude that the log-quadratic relation in Eq[T]is not the 
result of increased scatter, i.e. outliers, at the high-mass end of the M^-n relation (see Graham & Driver 2007a for more details). 




Figure 5. Logarithmic difference between the r- and g-band 
Sersic indices from individual galaxies in the SDSS-VAGC (Blan- 
ton et al. 2005), in which a single R 1 /™ function was fitted to 
each galaxy and the Sersic index restricted to lie between 0.2 and 
6. The peak of the distribution for galaxies with g-band Sersic 
indices greater than 2.0 (i.e., predominantly bulge-dominated, 
early-type galaxies) is 0.003 dex. The peak of the distribution for 
galaxies with n g < 2.0 (i.e., predominantly late-type galaxies) is 
0.06 dex. The choice of n g = 2.0 originates from the Sersic index 
bimodality plots in Driver et al. (2006a). While the right-hand 
panel is based on all the SDSS-VAGC data, the left hand panel 
shows a random sample of 75,000 with n g > 2.0 and 75,000 with 
n g < 2.0. Such a restriction prevents washing away the structure 
in the left panel with too many data points. 



-2 
-3 
-4 
-5 

0.5 


-0.5 

-2 

-3 
-4 
-5 

0.5 


-0.5 

-2 

-3 
-4 
-5 

0.5 


-0.5 




All types 



Early types 




Late types 



10 B 6 

log(M BH /M Sun ) 

Figure 6. Black hole mass functions for our three (colour) 
samples, e.g. S3 = Sample 3, see Section 15.21 Early-type galax- 
ies are those with B/T > 0.4 (see FiggJ. The best- fitting 3- 
parameter models (equation [2J are shown. Model parameters 
are given in Table \3\ The number density shown along the ver- 
tical axis is expressed in units of /i 3 Mpc — 3 per decade in 
SMBH mass. The fitting is performed over the SMBH mass range 
10 9 > M^/Mq > 10 6 . Residuals about the models are shown 
beneath each fit. 
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Figure 7. Cumulative SMBH space density, Pbh- Our early- type 
galaxies are those with B/T > 0.4, while our late-type galaxies 
have 0.01 < B/T < 0.4. The cutoff at an absolute magnitude 
of —18 B-mag marks the boundary to which we can trust our 
data. The horizontal shading is centred on the (Ho corrected, i.e. 
41 per cent increased) solution for (early+late)-type galaxies by 
Shankar et al. (2004, their section 3.2), and shows their quoted 
l-cr uncertainty. 
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Figure 8. Our observed SMBH mass functions, along with our 
fits to Sample 3 (solid curves), overlaid with others estimates 
of the mass function. Shown in the top panel arc curves from 
Shankar et al. (2004; S04) and Aller & Richstone (2002): (AR02)a 
is based on the Madgwick et al. (2002) luminosity function, and 
(AR02)b is based on the Marzke et al. (1994) luminosity function. 
Shown in the lower panel are curves from Shankar et al. (2004), 
and also McLure & Dunlop (2004; MD04) who used two inde- 
pendent techniques (see Section [1J. The (Mbh-o")-derived mass 
functions can be seen to turn down at low-masses, in agreement 
with our (Mbh-^)-derived mass function. We have corrected all 
curves to Ho = 70 km s — 1 Mpc -1 , following the prescription 
given by Graham & Driver (2007b). 
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Figure 9. Our (extreme) high-mass SMBH mass function con- 
structed using the linear M^^-n relation which excluded the three 
highest-n galaxies (see Section l5.1l l. The solid curves are the opti- 
mal fits shown in Figure [6] The dashed curve in the middle panel 
is the (-Mbh - L)-derived mass function given in McLure & Dunlop 
(2004). 
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Figure 10. The (u — r) c core colour distribution of the three 
spheroid samples. Only galaxies with B/T > 0.01 and Mg < —18 
mag are shown. A colour cut at (u — r) c = 2.00 mag reflects the 
central minimum observed in the bimodal distribution of core 
colours shown in Liskc ct al. (2007, in prep.). 



